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Beschreibung 

Verfahren und Anordnung zur Vorhersage von Mefidaben anhand 
vorgegebener Mefidaten 

Die Erfindung betrifft ein Verfahren und eine Anordnung zur 
Vorhersage von Mefidaten anhand vorgegebener Mefidaten. 

Ein technisches System wirft oftmals die Frage nach einer 
Prognose anhand bekannter (MeB-)Daten auf, insbesondere im 
Hinblick auf eine Fehleranf alligkeit oder eine 
Kostenabschatzung. 

Ubernimmt eine solche Prognose ein Experte, so ist diese 
Prognose zumeist f ehlerbehaf tet . Eine exakte Bestimmung kann 
der Experte, zumindest fur ein hinreichend komplexes System, 
nicht vornehmen. 

Ein stochastischer PunktprozeB, insbesondere ein Poisson- 
Prozeli, ist aus [1] bekannt. 

Die Aufgabe der Erfindung besteht darin, automatisch eine 
Vorhersage (Prognose) von Melidaten anhand vorgegebener 
Mefidaten zu ermoglichen. 

Diese Aufgabe wird gemali den Merkmalen der unabhangigen 
Patentanspruche gelost. Weiterbildungen der Erfindung ergeben 
sich auch aus den abhangigen Anspruchen, 

Zur Losung der Aufgabe wird ein Verfahren zur Vorhersage von 
Mefidaten anhand vorgegebener Melidaten angegeben, bei dem ein 
stochastischer Prozefi an die vorgegebenen Mefidaten angepafit 
wird. Ab einem vorgegebenen Zeitpunkt werden Simulationslauf e 
bis zu einem Endzeitpunkt durchgef uhrt . Fur jeden 
Simulationslauf werden die prognostizierten Melidaten 
bestimmt. Fur den Endzeitpunkt erfolgt die Vorhersage von 
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Mefidaten innerhalb eines Wertebereichs, der durch die 
prognostizierten Mefidaten bestimmt ist. 

Eine Weiterbildung besteht darin, dafi ein Konf idenzintervall 
fiir die Vorhersage von Mefidaten bestimmt wird, indem die a% 
kleinsten und die b% grofiten prognostizierten Mefidaten 
eliminiert werden. Insbesondere kann a% = b% sein. 
Beispielsweise kann somit ein 95%-iges Konf idenzintervall 
bestimmt werden, indem die 2,5% kleinsten und die 2,5% 
grofiten prognostizierten Mefidaten unberucksichtigt bleiben. 

Ein Vorteil besteht darin, dafi von einem vorgegebenen 
Zeitpunkt aus die Mefidaten mit einer in einem 
Konf idenzintervall liegenden Genauigkeit vorhergesagt 
(prognostiziert)- werden konnen. Dies ermoglicht bereits in 
einem fruhen Stadi um die Erkennung beispielsweise einer 
ErfullbaB^keit bzw. Nichterf ullbarkeit einer mit den Mefidaten 
verbundenen Aufgabe, so dafi daran geeignete Mafinahmen 
geknlipft werden konnen, um einer prognostizierten 
Nichterf ullung entgegenzuwirken . Dies ist insbesondere von 
Bedeutung bei einem komplexen System, z.B. einem 
Softwareentwicklungsprozefi, bei dem innerhalb einer spateren 
Testphase aufzeigbar ist, inwieweit eine Zeitplanung bis zur 
Fertigstellung der Software eingehalten werden kann. Umso 
wichtiger ist es, hierbei friihzeitig einer deutlichen M 
Verzogerung, beispielsweise in einer Integrationstest-Phase,* 
entgegenwirken zu konnen. Dies wirkt sich zum einen auf die 
Erfullbarkeit innerhalb der gesetzten Frist (Zeitrahmen) , zum 
anderen direkt auf die Kosten aus, da eine Nichterf ullung in 
der vereinbarten Zeit regelmafiig zusatzliche Kosten 
verursacht, 

Eine Ausgestaltung besteht darin,. dafi der stochastische 
Prozefi ein nichthomogener Poisson-Prozefi ist. 

Insbesondere ist es eine Ausgestaltung, dafi die Mefidaten 
Anzahlen von Fehlern sind. Dies entspricht beispielsweise 
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einer Sof twareentwicklung, bei der abhangig von in einer 
Testphase gemessenen Fehlern ein Reifegrad derselben 
dokumentiert wird. Abhangig von diesem Reifegrad resultiert 
direkt die Fertigstellung, d.h. solange nicht ein Grofiteil 
der Fehler aus der Software entfernt wurde, kann diese nicht 
an Kunden ausgeliefert warden. Dies driickt sich insbesondere 
durch Aufwand (beim Testen und Korrigieren der Fehler) und 
Kosten (fur die Zeitverzogerung bei der Auslief erung) aus. 

Zur Losung der Aufgabe wird auch ein Verfahren zur Vorhersage 
von Mefidaten anhand vorgegebener MefSdaten angegeben, bei dem 
ein stochastischer ProzeiJ an die vorgegebenen Meiidaten 
angepaftt wird. Es wird ein Intervall ermittelt, indem anhand 
des stochastischen Prozesses gewonnene 
Wahrscheinlichkeitswerte der Groiie nach urn einen 
Erwartungswert sortiert werden. Die Vorhersage von Meiidaten 
erfolgt durch Orientierung an dem Intervall, insbesondere an 
den Wahrscheinlichkeitswerten innerhalb des Intervalls. 

Fine Weiterbildung besteht darin, dali die anhand des 
stochastischen Prozesses gewonnenen Wahrscheinlichkeitswerte 
der Groiie nach symmetrisch um den Erwartungswert sortiert 
werden. Damit ist insbesondere gemeint, daii der grofite 
Wahrscheinlichkeitswert die Mitte des Intervalls, also den 
Erwartungswert kennzeichnet , wohingegen der nachst groiiere 
Wahrscheinlichkeitswert rechts- bzw. linksseitig des 
Erwartungswertes angeordnet wird. Der nachfolgend nachst 
hohere Wahrscheinlichkeitswert wird gegenuber auf der anderen 
Seite des Erwartungswertes symmetrisch angeordnet. 

Auf diesem Wege erhalt man analytisch (konstruktiv) ein 
Intervall, das wiederum durch seine Breite angibt, welche 
Wahrscheinlichkeitswerte fur die Vorhersage der Meiidaten eine 
Rolle spielen. 

Insbesondere ist es eine Ausgestaltung, dali die Breite des 
Intervalls bestimmt wird, indem diejenigen 
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Wahrscheinlichkeitswerte unberucksichtigt bleiben, die 
unterhalb einer vorgegebenen Schranke liegen, 

Dadurch ergibt sich ein Intervall (Konf idenzintervall) , 
welches durch die Schranke eine bestimmte Breite aufweist. 
Diese Breite entspricht der Sicherheit der Vorhersage von 
Meftdaten. 

Geht man davon aus, daJi der stochastische Prozeli ein 
nichthomogener Poisson-Prozeli ist, so wird insbesondere auf 
einer Zeitachse t durch den nichthomogenen Poisson-Prozefi 
eine Schrittweite bestimmt, die angibt, wann der nachste 
Fehler auftritt. Durch die Eigenschaft der 
Gedachtnislosigkeit des nichthomogenen Poisson-Prozesses wird 
von jedem zu ei^Fiem bes-^mmten Zeitepunkt auf getretenen Fehler 
"gedacht.rsid5slos"^?mach^ve^ ^iderV^nactis^t'en -Fehler 

kenn z e i ©h^enderi ^^idii^p^^ ge siixcht . 

Auch wird**^ur 4^ der^ AufgabeHiei»ne*^novirdnung zur 

Vo r he r s a^ge von^M^®da t"en'i*^*.anhand ^vox^g^'g-^teene r" MeiJ da t en 
angegeben, bei der eine' Prozessoreinheit vorgesehen ist, die 
derart eingerichtet ist, dali 

a) ein stochastischer Prozeli an die vorgegebenen Meiidaten 
anpafibar ist; 

b) ab einem vorgegebenen Zeitpunkt Simulationslauf e bis M 
2u einem Endzeitpunkt durchftihrbar sind; ^ 

c) fiir jeden Simulationslauf die prognostizierten 
MeBdaten ermittelbar sind; 

d) fiir den Endzeitpunkt die Vorhersage von MeBdaten 
innerhalb eines Wertebereichs , der durch die 
prognostizierten MeBdaten bestimmt ist, vorhersagbar 
ist . 

Ferner wird zur Losung der Aufgabe eine Anordnung zur 
Vorhersage von MeBdaten anhand vorgegebener MeBdaten 
angegeben, bei der eine Prozessoreinheit vorgesehen ist, die 
derart eingerichtet ist, daB 
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a) ein stochastischer Prozefi an die vorgegebenen MeJidate 
anpaftbar ist; 

b) ein Intervall ermittelbar ist, indem anhand des 
stochastischen Prozesses gewonnene 
Wahrscheinlichkeitswerte der GroBe nach um einen 
Erwartungswert sortiert warden; 

c) die Vorhersage von Mefidaten innerhalb der Grenzen des 
Intervalls erf olgt . 

Die Anordnungen sind insbesondere geeignet zur Durchfuhrung 
der erf indungsgemalien Verfahren oder der vorstehend 
erlauterten Weiterbildungen . 

Ausfuhrungsbeispiele der Erfindung werden nachfolgend anhand 
der Zeichnung dargestellt und erlautert. 

Es zeigen 

Fig.l eine Skizze, die eine akkumulierte Fehleranzahl uber 
einem Testzeitraum darstellt; 

Fig, 2 eine Skizze, die uberlagerte Konf idenzintervalle fur 
verschiedene ProzeBmodelle darstellt; 

Fig. 3 ein Flu/idiagramm mit Schritten eines Verfahrens zur 

Vorhersage von Meftdaten anhand vorgegebener Meiidaten; 

Fig. 4 ein weiteres Flu/idiagramm mit Schritten eines 
Verfahrens zur Vorhersage von MeJSdaten anhand 
vorgegebener Meiidaten; 

Fig. 5 eine Prozessoreinheit ; 

Um eine zu erwartende Fehleranzahl in einem technischen 
Prozeli, beispielsweise in einem Sof tware-Entwicklungsprozefi, 
prognostizieren zu konnen, werden nichthomogene Poisson- 



I 

GR 98 P-5929 



6 

Prozesse (NHPP) kalibriert, d.h. an MeJidaten, z.B. das 
Auftreten der Fehler uber die Zeit, angepalit . Mit 



{N(t)}t?^RH 



(1) 



10 



15 



20 



wird ein zu dem stochastischen Punktprozeli (nichthomogenen 
Poisson-Prozeli) assoziierter Zahlprozefi und mit einem 
Zeitpunkt ein Ende eines Testzeitraumes, also ein 
Zeitpunkt an dem die vorgegebenen Daten enden> bezeichnet. Es 
werden die stochastischen Prozesse 



{u(t)}t 
{o(t)}t 



und 



gesucht mit 



p(u(t)^ <^N(t) ^^o)^^(t)ifcN(£o) =^^no) ma 



(2) 



(3) 



(4) , 



fur all# *Zei%punk^^ t>to und gegebene Werte a e (o,l) 
(Konf idenzndveau) sowie ng e N . Nachfolgend werden 
insbesondere die Zuwachse des stochastischen Zahlprozesse 
bezogen auf den Zeitpunkt tg betrachtet. 



25 Fur den hier vorliegenden Fall, in dem Gleichung (1) einen 
nichthomogenen Poisson-Prozeli bezeichnet, gilt (vgl- [1]) 



M^l) - «(to) = ^) = exp(- [i(ti) - i(to)]) . 



(5) 



30 



f iir 



0 < to < ti < 00, e B Nq 



(6) 



35 



und eine Intensitat (englisch: intensity, mean measure, mean 
value function) 
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i: R"*" ^ R''",t »-» i(t>= EN(t) 



(7) 



Da gemali der Eigenschaft des Poisson-Prozesses die Zuwachse 
(hier: Fehlerzuwachse) unabhangig von vorangegangenen 
Zuwachsen sind, kann Gleichung (5) fiir Zeitpunkte t>to zur 
Bestimmung eines (moglichst kleinen) Intervalls 



[gu' go] = [gu(t), go(t)] c Nq 



(8) 



vereinfacht werden zu 



go 



2 p(N(t) - N(to) = ^) > a 

^=gu 



(9) 



Aufgrund der Unimodalitat der Poisson-Zahldichte kann ein 
Intervall [guf go] folgt ermittelt werden: 

Schritt 1: sortiere die Elementarwahrscheinlichkeiten 

: = p(N(t) - N(to) = e), e No 

in absteigender Reihenfolge und bezeichne die 
derart sortierten Werte itiit 



P(0)' P(l)' 



(d.h. {po/ PI/ . • .} = {p(0)' P(l)' • • •) 

P(o) ^ P(i) ^ . • . ) ; 



Schritt 2: bestimme £ 



min 



= min 



i e No 



E P(i) ^ « 

i = 0 



Schritt 3: bestimme eine Indexmenge 



I : = |io, . . . , ifjnin} ^ ^0 mit 
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} = K> ^(<min)}' 



Schritt 4: setze- : = min{i} und go : = max{i} . 

iel iGl 



Das Interyall aus Gleichung .(8) wi-rd^auch als Prognose- 
Intervall bezeichnet. 

S'bochastlsche Simulation (zwel'ber Ansa'bz) 

Das beschriebene Konf idenzintervall kann simulativ bestimmt 
werden durch folgende Schritte: 

Schritt 1: s4i:artei»arfr-^ed^t^^ 

mxm # N*^^©i^d*i^aMde3^^ 

^^r^fdti f^'^(!iempig^4^/l^aliMit^e rl^^B^ o^e'lkrv^de J. ]&^^b^s»fke r e n ; 

Schritt 2: ^%eetide??^eid«ne#^S^ 

gewlinseh^e-^l}n;d^eitpuo?i^^ erreicht ist; 

Schritt 3: wiederhole Schritt 2 solange bis alle 
Simulationslauf e beendet sind; 




Simulationslauf generierten Fehler in dem 
Zeitraum (to/teJf i=l,...,in, in absteigender 
Reihenfolge und bezeichne die derart sortierten 



Werte mit N(-L)(te), . . - , N(j^)(te) 



Schritt 5: 



setze 



gu • = 

go - = 



N(Lm.a/2j)(^e) und 
^(rm.(l-a/2)l)(^e). 
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d.h. eliminiere die (lOO • (l - a) / 2)% kleinsten 
und groISten Werte. 

5 Daraus resultiert unmittelbar das Konf idenzintervall . 



Jeder einzelne Simulationslauf basiert auf einem 
Simulationsalgorithmus, der aus (vgl. [2]) bekannt ist; 
Simulative Erzeugung von Zwischenankunf tszeiten eines 
nichthomogenen Poisson-Prozesses : 

Schritt 1: Setze X:= sup {A,(t)} , wobei gilt: 

t>ts 

. . di 

m : = - 



(lo; 



Schritt 2: erzeuge eine mit dem Parameter X exponential- 
verteilte (Pseudo-) Zuf allsvariable X, d.h. 
X: = -log(u)/A^, wobei U auf (0,1) gleichverteilt 

ist . 



0 Schritt 3; 



Schritt 4 : 



erzeuge eine auf (0,1) gleichverteilte 
Zuf allsvariable U. 

falls U < A,(ts x) / X , dann setze t* : = tg + X ; 
andernfalls setze tg : = tg + X und gehe zu 

Schritt 1. 



Fig.l zeigt beispielhaft eine Skizze, die eine akkumulierte 
Fehleranzahl tiber einem vorgegebenen Testzeitraum darstellt, 
Ab einem Zeitpunkt tg zeigt sich ein Vorhersageintervall fur 
alle Zeitpunkte tg + x. 



Allgemein ergibt sich die Ableitung der Intensitat i gemaft 
Gleichung (10) zu X. Beispielsweise ergibt sich: 
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a) kit) = a - b • c • exp(- bt*^) • t^""^ 

( A.(t) ist fiir c < 1 streng monoton f aliened, fiir c > 1 
unimodal mit einem eindeutigen Maximum an einer Stelle 



10 In Fig. 2 ist eine Skizze gezeigt, die uberlagerte 

Konf idenzintervalle darstellt. Insbesondere wird dadurch 
ersichtlich^ daB eine etwaige Prognose umso breiter streut^ 
je* weiter sie. in die Zukunft reicht. Insbesondere konnen 
durch verschiedene ProzeBmodelle berechnete 

15 Konf idenzintervalle auf. die in Fig . 2 dargestellte Art 
visualisiert^ werden. 

In Fig. 3 ist ein Flufidiagramm mit Schritten eines Verfahrens 
zur Vorhersage von MeBdaten anhand vorgegebener MeBdaten 

20 dargestellt. In einem Schritt 301 wird ein stochastischer 

ProzeB, insbesondere ein nichthomogener Poisson-Prozeli (als 
Vertreter fiir einen stochastischen Zahlprozeli) an vorgegebe 
Melidaten angepalit. In einem Schritt 302 werden 
Simulationslauf e durchgef uhrt , die ausgehend von dem 

25 Zeitpunkt tg bis zu einem zu prognostizierenden Endzeitpunkt 
te laufen. Fur jeden Simulationslauf werden in einem Schritt 
303 prognostizierte Melidaten bestimmt und eine Vorhersage von 
Melidaten auf einen Bereich eingeschrankt , der von den durch 
die Simulationslauf e bestimmten Meiidaten abgedeckt wird 

30 (siehe Schritt 304) . In einem Schritt 305 wird ein 

Konf idenzintervall bestimmt, indem jeweils ein vorgegebener 
Teil groiiter und kleinster prognostizierter Meiidaten 
unberiicksichtigt bleibt (dies entspricht besagtem Bereich) . 
In einem Schritt 306 wird das Verfahren beendet. 




5 



X ergibt sich nach den obigen Oberlegungen zu 
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Fig. 4 zeigt ein weiteres FluBdiagramm mit Schritten eines 
Verfahrens zur Vorhersage von MeBdaten anhand vorgegebener 
Meftdaten, Dazu wird in einem Schritt 401 eine Anpassung eines 
stochastischen Prozesses, insbesondere eines nichthomogenen 
Poisson-Prozesses an die vorgegebenen Mefidaten durchgef iihrt . 
Anhand des stochastischen Prozesses werden 
Wahrscheinlichkeitswerte bestimmt, die der GrofSe nach um 
einen Erwartungswert sortiert werden (vergleiche Schritt 
402) . Durch diesen Sortiervorgang wird ein Intervall, hier 
ein Konf idenzintervall, bestimmt. Die Breite des 
Konf idenzintervalls ergibt sich durch Vergleich der 
akkumulierten Wahrscheinlichkeiten mit einer vorgegebenen 
Schwelle. Durch das Konf idenzintervall ist, wie oben erklart 
wurde, eine Verteilung bzw. Unscharfe von einem Zeitpunkt tg 
in die Zukunft gegeben, die eine Abschatzung der Meiidaten in 
der Zukunft ermoglicht (siehe Schritt 403). In einem Schritt 
404 wird das Verfahren beendet. 

In Fig. 5 ist eine Prozessoreinheit PRZE dargestellt. Die 
Prozessoreinheit PRZE umfafit einen Prozessor CPU, einen 
Speicher SPE und eine Input /Output-Schnittstelle lOS, die 
iiber ein Interface IFC auf unterschiedliche Art und Weise 
genutzt wird: Uber eine Graf ikschnittstelle wird eine Ausgabe 
auf einem Monitor MON sichtbar und/oder auf einem Drucker PRT 
ausgegeben. Eine Eingabe erfolgt uber eine Maus MAS oder eine 
Tastatur TAST. Auch verfiigt die Prozessoreinheit PRZE uber 
einen Datenbus BUS, der die Verbindung von einem Speicher 
MEM, dem Prozessor CPU und der Input/Output-Schnittstelle lOS 
gewahrleistet . Weiterhin sind an den Datenbus BUS zusatzliche 
Komponenten anschlielibar , z.B. zusatzlicher Speicher, 
Datenspeicher (Festplatte) oder Scanner. 

Nachfolgend werden ein Algorithmus zur Bestimmung von 
Konf idenzintervallen fur Prognosen und ein Algorithmus zur 
simulativen Bestimmung von Konf idenzintervallen fur Prognosen 
in der Notation der Programmiersprache C angegeben. 
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Programm 1 : 

/* Bestimmung von Konfidenzintervallen fQr Prognosen */ 
/* basierend auf dem verallgemeinerten Goel-Okomoto-Modell */ 
5 #include <stdlib.h> 
#inelude <math.h> 
#include <stdio.h> 

#define true 1 
10 #define. false -1 

double mv _jgenGO(double,double,double,double); double poisson(double,long); void ki_nhpp(); 
int main(argc,argv) 
Int argc; 
1 5 char *argv[]; 
{ 

double a,b,c,bt.st,kn; 
long low.uppjauf; 

20 if(argc<7){ 

printf("\n\nZuwenig Argumente! \n\n"); 

printf('7Kufruf:.^s <Par1> <Par2> <Par3>^Startzeit> <Endzeit>", 
"<KNiveau>\n\n"- argv[0])f-return 1 ; 
} 

25 

a = at0f(arg^O1]); 
b = atof(arg>/[2]); 
c = atof(argv[3]); 
bt= atof(argv[4]); 
30 st= atof(aFgV[5]); 
kn= atof(argv[6]); 

for (lauf=1;lauf< ;lauf++) { 

ki_nhpp(mv_gehGO,a,b.c,bt,bt+lauf*(st-bt)/10.,kn,&low,&upp); 
35 printffZeitpunkt: %8.2f Fehlerintervall: [%d,%d]\n", 
bt+laur(st-bt)/10., low. upp); 
} 

return 0; 
} 

40 

double mv_genGO(x,a,b,c) 
double x,a.b,c; 

{ return( a*(1.0-exp(-b*pow(x,c))) ); } 

4 5 double-(poisson(lambda,wert) 
double lambda; 
long wert; 
{ 

long i; 
50 double itval.hv; 

if (lambda<600) { 
itval = exp(-lambda); 

for (i=wert;i>=1;i") { itval *= lambda/(double)i; } 
55 } 
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else { 

hv = xp(-lambda/(doubIe)wert); 
itval=1.0: 

for (i=wert;i>=1;i-) { itval *= lannbda/(double)i*hv: } 

5 } 

return ( itval ); 

} 

void ki_nhpp(mv_nhpp. par1_nhpp, par2_nhpp, par3_nhpp, 
10 start^time, stop_time, k_niveau, lower, upper) double mv_nhpp(double,double,double.double); 

double par1_nhpp, par2_nhpp, par3_nhpp. start_time. stop_time. k_niveau; long *lower, *upper; 

{ 

long lauf; 

int lborder,modJow,mod_upp; 
1 5 double sum,tnnp_mv, vaM, va!_u; 

tmp_mv = mv_nhpp{stop_time.par1_nhpp,par2_nhpp,par3_nhpp) - 
mv_nhpp(start_time,par1_nhpp,par2_nhpp.par3_nhpp); lauf = (long)tmp_mv; 
*lower = lauf-1 ; 

*upper = lauf+1; modjow= false; mod_upp= false; sum = poisson(tmp_mv,lauf); valj = 

poisson(tmp_mv,*lower); vaLu = polsson(tmp_mv,*upper); 
while (sum<k_niveau) { 
if (valj<val_u) { 
sum += vai_u; 
25 (*upper)++; 

Iborder = false; 
mod_upp = true; 

val_u = poisson(tmp_mv,*upper); 
} 

30 else{ 

sum += valJ; 
(*lower)-; 
Iborder = true; 
modjow = true; 
35 valj = poisson(tmp_mv,*iower); 

^^^V if (Iborder == true) { (*lower)++; } 

^^Wl) else { (*upper)-; } 

if (modjow == false) { (*lower)++; } 
if (mod_upp == false) { (*upper)-; } 

4 5 retum; 
} 
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Programm 2 : 

/* Simulative Bestimmung von Konfidenzintervallen fOr Prognosen */ 
/* basierend auf dem verallgemeinerten Goel-Okomoto-Mode!l */ 
5 #include <stdlib.h> 
#include <math.h> 
#include <time.h> 
#include <stdio.h> 
#include <values.h> 

10 

#define true 1 
#define false -1 

double drand48(void); 
15 void srand48(long); 

double sim_exp(double); double lambda_genGO(double,double,double,double); void sim_nhpp(); 
int main{argc,argv) 
int argc; 
2 0 char *argv[]; 
{ 

timej t; double a,b^e;bt,stvpnt[1000000];check^time[12]; long lauf,no_pnt,seed_run; int clauf; 
FILE *datei; 
if (argc<6) { 

2 5 printf("\n\n2uwenig Argumente! \n\n"); 

printfC7\ufruf:^s^^<Par1 > <Par2> <Par3> <Start2eit> <Endzeit>\n\n", 
argv[0]); return 1; 
} 

3 0 datei = fopen{"sinn:seed"/r); 

if (datei==NULL) { 
seed_run = 1 ; 
} 

else{ 

3 5 fscanf(datei,"%6d",&seed_run); fclose(datei); seed_run++; 



datei = fopen("sim.seed"."w+"); 
fprintf(datei. "%6d\n", seed_run ); 
4 0 fclose(datei); 

time (&t) ; /* Initialisierung des 7 

t +T= seed_run*100 ; /* Zufallszahlengenerators */ 

srand48 ((unsigned long) t) ; /* mit Hilfe der Systemzeit V 
4 5 a- atof(argV[1]); 

b = atof(argv[2]); 

c = atof(argv[3]); 

bt= atof(argv[4]); 

st= atof(argv[5]); 

50 

sim_nhpp(lambda_genGO.a,b.c,bt,st.&pnt.&no_jDnt); 
for (lauf=1;lauf<=no _pnt;lauf++) { 
printf("%15.7f %10d \n", pnt[lauf]. lauO; 

55 
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datei = fopen("ki.tmp","a"); 
for (lauf=1;lauf< ;lauf++) { 
check_time[lauf] = bt+laur(st-bt)/10.; 

} 

check Jime[1 1] = pnt[noj)ntl+1 ; T grdBer als die grdQte 
simulierte Zeit */ 
clauf = 1 ; 

for (lauf=1;lauf<=no_j)nt;lauf++) { 

while (pnt[lauf]>=check_time[claufl) { fprintf(datei. "%8.2f %6d check_time[clauf], lauM); clauf++: 
} 

} 

if {pnt[noj)nt] <check_time[10]) { 

for {lauf=clauf;lauf< ;lauf++) { fprintf(datei, "%8.2f %6d check Jime[lauf]. nojtnt); 
} 

} 

fprintf(datei. 'An"); 
fclose(datei); 

return 0; 
} 

double sim_exp(lambda) 
double lambda; 

{ return( -log(drand48())/lambda ); } 

double lambda_genGO(x,a,b.c) 
double x,a,b,c; 

{ return( a*b*c*pow(x,c-1)*exp(-b*pow(x,c)) ); } 

void sim_nhpp(lambda_nhpp. par1_nhpp, par2_nhpp, par3_nhpp, 

startjime, stop_time, path, no_points) double lannbda_nhpp(double,double,double,double); double 
par1_nhpp, par2_nhpp, par3_nhpp. start_time. stop_tinne; double path[]; long *no_points; 

{ 

double sim_time,x,u,x_bar,lambda_bar; 
*no _points=0; 
sim_time = start_time; 

do{ 

if {par3_nhpp<=1) { lambda_bar = Iambda_nhpp(sim_tlme,par1_nhpp.par2_nhpp,par3_nhpp); 
} 

else { 

x_bar = pow((par3_nhpp-1.0)/par2_nhpp/par3_nhpp,1.0/par3_nhpp); 
if {sim_time>=x_bar) { 

lambda_bar = Iambda_nhpp(sim_time,par1_nhpp,par2_nhpp,par3_nhpp); 
} 

else { 

lambda^bar = Iambda_nhpp(x_bar,par1_nhpp,par2_nhpp,par3_nhpp); 
} 

} 

X = sim_exp(lambda_bar); 
u = drand48(); 
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if (u<=lambda_nhpp(sim_time+x,par1_nhpp.par2_nhpp,par3_nhpp)/lambda_bar) { (*no_points)++; 
path[*noj30ints]=sim_time+x; 

} 

sim_time+=x; 
5 } 

while (sim_time'^:=stop_time); 

return; 

} 
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Progranm 3 : 

r Bestimmung der Konfidenzintervalle aus den Simulationsdaten */ 
r (die Simulationsdaten werden in aufsteigender Reihenfolge sortiert) */ 
5 #include <stdlib.h> 
#include <math.h> 
#include <stdio.h> 

int qsortJcmp(int*,int*); 
1 0 int qsortjcnnp(x.y) 
int *x. *y; 
{ 

if (*x<*y) { return ( -1 ); } 
else if rx==*y){ return ( 0);} 
15 else { return ( 1 ); } 

} 

int main(argc,argv) 
int argc; 
2 0 char*argv[]; 

{ 

int pnt[11][1 00000]; 
int qs[1 00000]; 
char *dname; 
25 intfrac.i; 

long lauf,lower_bound,upper_bound; 
long l,no_pnt,seed_run; 
double ctime[11],x; 
FILE *datei; 

30 

if(argc<3){ 

printfCAnXnZuwenig Argumente! \n\n"); printf("Aufruf. %s <Dateiname> <Konfidenzniveau (ii 
%%)>\n\n". argv[0]); return 1; 

35 } 

dname = argv[1]; 
frac = 100-atoi(argv[2]); 
lauf = 0; 

40 

datei = fopen(dnanne,V); 
if (datei==NULL) { return 1 ; } 
else { 

while (!feof(datei)) { 
4 5 lauf++; 

for(i=1;i<=9;i++){ 

fscanf(datei;%8lf %6d". &ctime[il, &pnt[i][lauf]); 

fscanf(datei."%8lf %6d \n". &ctime[10]. &pntt10][lauf]); 
50 } 

fclose(datei); 

} 

lower_bound = (long)floor(lauf*frac/200.); upper_bound = (iong)ceil(lauf*(200.-frac)/200.); 
55 if (lower_bound<1) {lower_bound = 1;} 
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printf("\n\n%2d%%-Sicherhettsbereich bei %d Simulationslaufen\n\n" 
100-frac,lauf); 
for (i=1;i< ;i++) { 
for (l=1;l<=lauf;l++){ 
5 qs[l] = pnt[i][l]: 
} 

qsort(&qs[1], lauf, sizeof(int), Aqsortjemp); 
printfC'Zeitpunkt: %8.2f Fehlerintervall: [%d,%cl]\n". 
0 ctime[i], qs[lower_bound], qs[upper_bound]); 



15 



return 0 ; 
} 
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Patentanspruche 

1. Verfahren zur Vorhersage von Meiidaten bis zu einem 
Endzeitpunkt anhand vorgegebener-Mefidaten, 

a) bei dem ein stochastischer Prozeli an die vorgegebene 
Melidaten angepaBt wird; 

b) bei dem ab einem vorgegebenen Zeitpunkt 
Simulationslaufe des stochastischen Prozesses bis zu 
dem Endzeitpunkt durchgefUhrt werden; 

c) bei dem fiir jeden Simulationslauf die prognostiziert. 
Mefidaten bestimmt werden; 

d) bei dem die Vorhersage von Melidaten durch Angabe ein 
Wertebereichs, der durch die prognostizierten MeBdate 
bestimmt ist, erfolgt. 

Verfahren nachv^Anspruch 1, 

bei dem ein. Konfidenzintervall fiir die Vorhersage von 
MeBdaten bestimmt , wird, indem die a% kleinsten und die b- 
grofiten prognostizierten Melidaten eliminiert werden. 

Verf ahren ^nach Anspruch 2, 
bei dem a% gleich b% ist. 

Verfahren nach einem der vorhergehenden Anspruche, 
bei dem der stochastische ProzeJi ein nichthomogener 
Poisson-Prozeli ist. 

Verfahren nach einem der vorhergehenden Anspruche, 
bei dem die Melidaten Anzahlen von Fehlern sind. 



6. 



Verfahren zur Vorhersage von Melidaten anhand vorgegebener 
Melidaten, 

a) bei dem ein stochastischer Prozeli an die vorgegebenen 
Melidaten angepalit wird; 

b) bei dem ein Intervall ermittelt wird, indem anhand des 
stochastischen Prozesses gewonnene 



98 P 5929 



21 

Wahrscheinlichkeitswerte der Grolie nach urn einen 
Erwartungswert sortiert werden; 
c) bei dem die Vorhersage von Melidaten innerhalb der 
Grenzen des Intervalls erfolgt. 

Verfahren nach Anspruch 6, 

bei dem die anhand des stochastischen Prozesses 
gewonnenen Wahrscheinlichkeitswerte der GroBe nach 
symmetrisch urn den Erwartungswert sortiert werden. 



Anordnung zur Vorhersage von Melidaten bis zu einem 
Endzeitpunkt anhand vorgegebener Melidaten, 
bei der eine Prozessoreinheit vorgesehen ist, die derart 
eingerichtet ist, daJi 

a) ein stochastischer Prozeli an die vorgegebenen Melidaten 
anpalibar ist; 

b) ab einem vorgegebenen Zeitpunkt Simulationslauf e des 
stochastischen Prozesses bis zu dem Endzeitpunkt 
durchfiihrbar sind; 

c) fur jeden Simulationslauf die prognostizierten 
Melidaten ermittelbar sind; 

d) die Vorhersage von Melidaten durch Angabe eines 
Wertebereichs, der durch die prognostizierten Melidaten 
bestimmt ist, erfolgt. 



Anordnung zur Vorhersage von Melidaten anhand vorgegebener 
Melidaten, 

bei der eine Prozessoreinheit vorgesehen ist, die derart 
eingerichtet ist, dali 

a) ein stochastischer Prozeli an die vorgegebenen Melidaten 
anpalibar ist; 

b) ein Intervall ermittelbar ist, indem anhand des 
stochastischen Prozesses gewonnene 
Wahrscheinlichkeitswerte der Grolie nach um einen 
Erwartungswert sortiert werden; 



■Ik 
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c) die Vorhersage von Mefidaten innerhalb der Grenzen de 
Intervalls erfolgt. 
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Zusammenf assung 

Verfahren und Anordnung zur Vorhersage von Mefidaten anhand 
vorgegebener Mefidaten 

Es werden ein Verfahren und eine Anordnung zur Vorhersage von 
Mefidaten anhand vorgegebener Mefidaten angegeben, bei dem/bei 
der ein stochastischer Prozefi an die vorgegebenen Mefidaten 
angepafit wird. Ab einem vorgegebenen Zeitpunkt werden 
Simulationslauf e bis zu einem Endzeitpunkt durchgef uhrt . Fur 
jeden Simulationslauf werden die prognostizierten Mefidaten 
bestimmt. Fur den Endzeitpunkt erfolgt die Vorhersage von 
Mefidaten innerhalb eines Wertebereichs, der durch die 
prognostizierten Mefidaten bestimmt ist. 
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